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Abstract 

We present a preliminary measurement of the branching fraction of the B meson decay 
B^ a5*"(1260)7r- with a5*"(1260) Tr+vr+vr". The data were recorded with the BABAR detector 
at the SLAC B factory PEPTI and correspond to 124 x 10*^ BB pairs produced in e^e an- 
nihilation through the T(4S') resonance. We find the branching fraction B{B^ a^(1260) 7r~) = 
(42.6 ±4.2 ±4.1) X 10"^. The fitted values of the ai(1260) parameters are niai = 1-19 ^ 0-02 
GeV/c2 and = 312 ± 55 MeV/c^. 
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1 INTRODUCTION 



We report on the preliminary measurement of the branching fraction a^(1260)7r~ with 

a^(1260) — > 7r"*'7r^7r~ T . The ai(1260) 3tt decay proceeds mainly through the intermediate 
states (7r7r)p7r and (7r7r)o-'/r. 

The study of this decay mode is complicated by open questions on the parameters of the ai (1260) 
meson. There are large discrepancies between these parameters when comparing results from 
analyses involving hadronic interactions [Ij and r decays jSj. Therefore, it is important to verify 
the theoretical prediction of the branching fraction for this decay mode and have new measurements 
of the ai(1260) parameters. It is also important to note that the a^(1260)7r^ channel can 

be used to measure the Cabibbo-Kobayashi-Maskawa angle a of the Unitarity triangle 

There has been no experimental observation of this decay mode. An upper limit of 49 x 10"^ at 
the 90% C.L. has been set by CLEO 5 for the branching fraction of B^ — > a^(1260)7r~, while the 
DELPHI IH] collaboration has set the 90% C.L. upper limit of 28 x 10~^ for the branching fraction 
of 5° ^ 4tt. 

Below we present the details of the analysis for the measurement of the branching fraction 
for B^ a^(1260)7r~ Att. Presently, we do not distinguish between the final states (7r7r)p7r 
and (7r7r)o-7r. Such an analysis would require a study of the angular distributions of the decay 
products. Background contributions from B^ decays to a2(1320)7r and 7r(1300)7r were assumed to 
be negligible. 

2 THE BaBar detector AND DATASET 

The results presented in this paper are based on data collected in 1999-2003 with the BABAR 
detector [7] at the PEP-II asymmetric e^e~ collider 8 located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 112 fb~^, corresponding to 124 million BB pairs, was recorded 
at the T(4S') resonance ( "on-resonance" , center-of-mass energy ^/s = 10.58 GeV). An additional 
12 fb~^ were taken about 40 MeV below this energy ("off-resonance") for the study of continuum 
background in which a light or charm quark pair is produced instead of an T(4S'). 

The asymmetric beam configuration in the laboratory frame provides a boost of /?7 = 0.56 to 
the T(45). Charged particles are detected and their momenta measured by the combination of a 
silicon vertex tracker (SVT), consisting of five layers of double-sided silicon detectors, and a 40- 
layer central drift chamber, both operating in the 1.5-T magnetic field of a solenoid. The tracking 
system covers 92% of the solid angle in the CM frame. 

Charged-particle identification (PID) is provided by the average energy loss {(SE/dx ) in the 
tracking devices and by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering 
the central region. A K/tt separation of better than four standard deviations [a) is achieved for 
momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta in the B decay final states. 
Photons and electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The EMC 
provides good energy and angular resolutions for detection of photons in the range from 30 MeV to 
4 GeV. The energy and angular resolutions are 3% and 4 mrad , respectively , for a 1 GeV photon. 

The flux return for the solenoid is composed of multiple layers of iron and resistive plate cham- 
bers for the identification of muons and long-lived neutral hadrons. 
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3 ANALYSIS METHOD 



Monte Carlo (MC) simulations jUI of the signal decay modes and of continuum and BB backgrounds 
are used to establish the event selection criteria. We select a||^(1260) candidates with the following 
requirement on the invariant mass in GeV/cP : 0.6 < ''^a^(i260) ^ ^'^ ' "^^^ intermediate dipion 
state is reconstructed with an invariant mass between 0.46 and 1.1 GeV/c^. 

We make several particle identification requirements to ensure the identity of the signal pions. 
For the bachelor charged track we require an associated DIRC Cherenkov angle between —2 a and 
+5 a from the expected value for a pion. 

A B meson candidate is characterized kinematically by the energy-substituted mass tties = 
^ {^s + Po • /Eq — and energy difference AE' = E*^ — \^/s, where the subscripts and B 
refer to the initial T(4S') and to the B candidate in the lab-frame, respectively, and the asterisk 
denotes the T(45) frame. We require |A^| < 0.2 GeV and 5.25 < mss < 5.29 GeV/c^. The 
momentum of a^(1260) in the center-of-mass frame is required to be between 2.3 and 2.7 GeV/c. 
To reduce fake B meson candidates we require p(x^) > 0.01 for the B vertex fit. The angular 
variable 7ia^ (cosine of the angle between the direction of the vr meson with respect to the flight 
direction of the B in the ai(1260) meson rest frame) is required to be between -0.85 and 0.85 to 
suppress combinatorics. 

To reject continuum background, we make use of the angle 9t between the thrust axis of the 
B candidate and that of the rest of the tracks and neutral clusters in the event, calculated in 
the center-of-mass frame. The distribution of cos 9t is sharply peaked near ±1 for combinations 
drawn from jet-like qq pairs and is nearly uniform for the isotropic B meson decays; we require 
I cos I < 0.65. The remaining continuum background is modeled from sideband data. We use 
Monte Carlo simulations of B^B^ and B^B^ decays to look for BB backgrounds, which can 
come from both charmless and charm decays. We find that the decay mode B^ — > D^n^, with 
D~ — > K^TT^TT^ and D~ K^7r~ , is the only significant background, and is included in the 
maximum likelihood fit. 

We use an unbinned, multivariate maximum-likelihood fit to extract the signal yield for B^ 
a^(1260)7r~. The likelihood function incorporates four uncorrelated variables. We describe the B 
decay kinematics with two variables: AE and rriEs, as mentioned above. We also include rua^ and 
a Fisher discriminant which describes the energy flow in the event. The Fisher discriminant 
combines four variables: the angles with respect to the beam axis, in the T(45) frame, of the B 
momentum and B thrust axis, and the zeroth and second angular moments -Lo,2 of the energy flow 
around the B thrust axis. The moments are defined by 

Lj = ^piX |cos6'i|-' , (1) 
i 

where 9i is the angle with respect to the B thrust axis of track or neutral cluster i, pi is its 
momentum, and the sum excludes tracks and clusters used to build the B candidate. 

Since the correlations between the observables in the selected data are small, we take the 
probability density function (PDF) for each event to be a product of the PDFs for the separate 
observables. The product PDF for event i and hypothesis j, where j can be signal, continuum 
background or BB background, is given by 

= VjiruEs) ■ V.iAE) . V,{:F) ■ V,{ma,). (2) 
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There is the possibiUty that a track from a signal event is exchanged with a track from the rest 
of the event. We call these events "self-cross-feed" (SCF) events. The fraction of SCF events with 
respect to the total number of signal events, fscF, is found to be 0.31 from Monte Carlo studies. 
The likelihood function for the event i is defined as : 

C = Usigil - fsCF)Vlig + nsigfsCpVscF + nqqVlg + UbbiPbEI + ^832^332 ' (3) 

where Ugig is the number of signal events, Uqq the number of continuum background events, J^g^i the 
number of BB background events D^ir^ with K^ir^ir^ and ngg2 the number of BB background 
events D^n^ with K'^tt^ . The extended likelihood function for all events is : 

where nj is the yield of events of hypothesis j found by the fitter, and N is the number of events 
in the sample. The first factor takes into account the Poisson fluctuations in the total number of 
events. 

We determine the PDFs for signal and BB backgrounds from MC distributions in each observ- 
able. For the continuum background we establish the functional forms and initial parameter values 
of the PDFs with data from sidebands in rriES or AE. We allow the signal rria^ PDF parameters 
and several background PDF parameters to float in the flnal fit. 

The distribution of rua^ in signal events is parameterized as a relativistic Breit-Wigner. The 
mES and AE distributions for signal are parameterized as double Gaussian functions. Slowly 
varying distributions are parameterized by linear functions. The combinatoric background in ttt-es 
is described by a phase-space- motivated empirical function ^j. We model the ^ distribution using 
a Gaussian function with different widths above and below the mean. 

Possible differences between Monte Carlo simulation and on-resonance data are investigated 
using the control sample B — > ■k~D^,D^ — > K~7r~^n^, which has a similar topology to the signal 
mode. 

4 SYSTEMATIC STUDIES 

Most of the systematic errors on the yields that arise from uncertainties in the values of the PDF 
parameters have already been incorporated into the overall statistical error, since they are floated 
in the fit. We determine the sensitivity to the other parameters of the signal PDF components by 
varying these within their uncertainties. The results are shown in the first row of Table ^ This 
is the only systematic error on the fit yield; the other systematics apply to either the efficiency or 
the number of BB pairs in the data sample. 

The uncertainty in our knowledge of the efficiency is found to be 0.8Nt%, where Nt is the number 
of signal tracks. We estimate the uncertainty in the number of BB pairs to be 1.1%. The fitting 
algorithm introduces a systematic bias of 3.1%, which was found from fits to simulated samples with 
varying background populations. Published world averages provide the B daughter branching 
fraction uncertainties. The systematic error from ai(1260)iir cross-feed background is estimated 
to be 5%, while the systematic error due to SCF is found to be 3%. A systematic error of 1% is 
assigned to potential contributions from i? — > 47r and B — > pmr. Finally, we account for systematic 
effects in cos^t (1%) and in the PID requirement (0.5%) on the prompt charged track. The values 
for each of these contributions are given in Tabled 
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Table 1: Estimates of the systematic errors (in percent). 



Quantity 


1 


Fit yield 


6.3 


Fit eff/bias 


1.7 


Track multiplicity 


1.0 


Tracking eff/qual 


3.2 


Number BB 


1.1 


SCF 


3 


aiK cross-feed 


5 


B — > 47r,p7r7r 


1.0 


MC statistics 


1.0 


cos 


1.0 


Total 


9.6 



5 RESULTS 

By generating (from the PDFs) and fitting simulated samples of signal and background, we verify 
that our fitting procedure is working properly. We find that the minimum ln£ value for the 
on-resonance data lies well within the ln>C distribution from these simulated samples. 

The efficiency is obtained from the fraction of signal MC events passing the selection criteria, 
adjusted for any bias in the likelihood fit. This bias is determined from fits to simulated samples, 
each equal in size to the data and containing a known number of signal MC events combined with 
events generated from the background PDFs. We find a fit bias of 0.97. 



Table 2: Final fit results. 



Fit quantity 


a5^(1260)7r- 


Fit sample size 




On-resonance 


32500 


Off-resonance 


2680 


Signal yield 




On-res data 


472 


Off-res data 
Selection e (%) 


c 9+10.8 
0-^-8.4 

19.4 


Track corr. 


0.953 


Fit-bias 


0.967 


US^ (%) 


50 


Stat. sign, (a) 


13.8 


e(xl0-6) 


42.6 ±4.2 ±4.1 







In Table |21 we show the results of the fits for on- and off-resonance data. We also show the 
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fitted signal yield, the efficiency (e), the daughter branching fraction product (H'^i)) the statistical 
significance, and the central value of the branching fraction. The statistical error on the number 
of events is taken to be the change in the central value when the quantity — 21n£ changes by one 
unit. The statistical significance is taken as the square root of the difference between the value of 
— 21n£ for zero signal and the value at its minimum. 

In Fig. n we show the mES; rua^, and rup projections made by selecting events with a 

signal likelihood (computed without the variable shown in the figure) exceeding a threshold that 
optimizes the expected sensitivity. 





Figure 1: Projections of mES (&)> (b), nia^ (c), and mass (d). Points with errors represent 
data, dotted lines the continuum and BB backgrounds, solid curves the full fit functions, and 
dashed curves the signal. These plots are made with a cut on the signal likelihood and thus do not 
show all events in the data sample. 

The fitted values of the ai(1260) parameters are = 1.19±0.02 GeV/c^ and Ta^ = 312±55 
MeV/c2 . 
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6 SUMMARY 



We have obtained a preliminary measurement of the branching fraction for meson decays to 
oj'"(1260) 7r~ with ai(1260) —>■ Stt. The measured branching fraction is: 

^(5° ^ a^(1260) TT-) = (42.6 ± 4.2 ± 4.1) x IQ-^, (5) 

where the first uncertainty is statistical and the second uncertainty is systematic. The fitted 
values of the ai(1260) parameters are niai = l-^^ ^ O-O^ GeV/c^ and = 312 ± 55 MeV/c^. 
These values are closer to those found in hadronic production of ai(1260) meson. 
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